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SUMMARY 


The  long-period  perturbations  of  the  orbital  elements  cf 
satellites  caused  by  direct  solar  radiation  pressure  are  considered 
for  satellites  of  high  "area-to-mass"  ratio.  A  series  of  orbits  is 
selected  to  illustrate  the  main  perturbation  effects.  The  graphs  are 
drawn  up  neglecting  planet  asphericity  and  the  effects  of  shadowing. 
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LONG-PERIOD  SATELLITE  PERTURBATIONS  _ 
UNDER  THE  ACTION  OF  SOLAR  RADIATION  PRESSURE  ' 


by 

E.  N„  Polyakhova 


Let  us  onsider  an  artificial  planetary  satellite  whose  orbital  position 
is  fixed  by  the  sic  Kepler  elements  (j  *  1,2,..., 6).  The  equation  of  the  per¬ 
turbed  motion  affected  by  light  radiation  pressure  can  be  written  as  an  equation 
in  a  small  parameter 


(tXl; 

____ 


'•a. 


K) 


/*=!.  2. 
\y  =  1.  2. 


(1) 


where  a,  the  ratio  of  the  modulus  of  the  perturbing  acceleration  to  the  gravita¬ 
tional  parameter  of  the  planet,  serves  as  the  small  parameter.  Equation  (1)  indi¬ 
cates  that  the  rijht-hand  terms  contain  the  orbital  elements,  the  sun’s  longitude 
and  the  satellite's  orbital  eccentricity  E. 

Equation  (1)  may  be  integrated  approximately  using  either  Picard's  method 
or  the  small -parameter  method,  since  for  an  equation  such  as  (1)  the  first  and  _ 
second  approximations  agree  for  both  methods  [lj.  The  transformation 

r  E 

J  9k  (4  >•©,  E)  dt  =  j'  -t..  (4  X3,  £)  dE,  (2) 

i.  ht 


where  x^  *  xj  (t  * 
ting  the  integral. 


tg)  are  the  initial  values  of  the  elements,  is  useful  in  evalua- 
Explicit  expressions  for  the  right-hand  terms  of  (1)  are  easily 


^  Translated  from  Vestnik  Leningradskogo  Universiteta  [Bulletin,  Leningrad 
University],  No.  7,  pp.  144-152  (1970). 
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obtained  from  the  known  differential  equation*  of  perturbed  motion,  the  equations 
being  in  Newtonian  form— -i.e.  containing  component*  of  the  perturbing  acceleration 
[2].  The  perturbing  acceleration  will  be  defined  by  the  formula* 


m 


F—  —  Fua,  F=r-Pr  4,  (3) 

where  i*  the  magnitude  of  the  solar  radiation  pressure  at  the  planet's  orbit, 

A/m  is  the  "surface-to-aass"  ratio  of  the  satellite,  u^  is  a  unit  vector  in  the 
earth*sun  direction.  Formulas  (3)  correspond  to  a  spherical  satellite  model,  for 
which  the  force  of  the  pressure  is  known  not  to  depend  on  the  reflectivity.  It  i3 
convenient  to  define  the  components  of  the  unit  vector  u®  in  an  orthogonal  system 
of  axes  located  at  the  pericenter  of  the  satellite's  orbit  (radial,  transverse  and 
bi-normai  directions).  The  derivations  of  these  formulas  may  be  found  in  reference 
C 3 J .  We  cite  here  the  final  expressions  for  the  projections  of  the  unit  vector 
on  these  axes: 


Formulas  for  the  perturbations  over  n  single  revolution  can  be  derived 
by  integrating  equation  (1),  taking  (2)  and  (4)  into  account,  with  the  condition 
that  the  orbital  elements  and  the  solar  longitude  can  be  considered  to  be  constants 
to  within  acceptable  accuracy  over  one  revolution  of  the  satellite: 


Ze=a'c(SQ,-TQt), 

U  —  —  a?s  W (Q*  cos  u.  —  Q,  sin  u>), 

22  =  —  —y  W  (Qi  cos  m  +  Qa  sin  u>). 


(5) 


(SQi  —  TQi)— cos  i2Q, 

ZAt0=?xa*seS  —  Vl~e:  (2«»  +  cos  i2Q)  -j- 

+  2a*3  Vr^F{SQi+TQl), 


where  the  quantities  Q^,  are  defined  by 


Qt~cos  £~  j cos 2c, 

Ql==??lW [3^“20  + e’Jsinff-f-JsinSf], 

Qa=  i(i_ei)C0S  2E, 

Qt=— 4-^1  —  ^(SE— 4<?$in£+y$ln2f?j, 
Qj=  -  7  ^3£ —  2*  sin  £  -  ±  sin  2£j , 

Qs  *=  —  e  cos  £  -f  ~  cos  2/.'. 


(6> 


Formulas  (5)  are  convenient  for  calculating  the  perturbations  over  any  number  of 
orbital  periods  q,  by  recalculating  the  elements  from  revolution  to  revolution  in 
accordance  with  the  scheme 

The  quantities  (6)  oust  be  calculated  within  the  limits  to E^,  which  means  that 
the  integration  is  carried  out  over  a  complete  period  of  change  in  E  from  0  to  2n, 
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froa  Which  we  exclude  the  region  of  shadowing  by  the  planet,  when  there  is  no 
acceleration  perturbation.  If  the  shadowing  limits  in  the  orbit  are  defined  as 
points  with  the  eccentric  anomalies  E^  (entering)  and  E2  (leaving),  the  integration 
must  be  carried  out  within  the  limits  0  to  F^,  and  froa  E2  to  2n  if  the  perigee 
(E  *»  0)  lies  in  light  (E2  -  >  0)  and  from  E2  to  if  the  perigee  lies  in  the 

shadow  (E2  -  E^  <  0). 

Formulas  (5)  for  the  case  of  total  illumination  can  be  used  for  calcu¬ 
lating  the  perturbations  over  one  or  a  number  of  revolutions,  or  also  over  some 
portion  of  a  revolution.  Having  obtained  these,  it  is  simple  to  follow  the  evo¬ 
lution  of  the  orbit  of  satellites  over  long  intervals  of  time  by  computing  the 
elements  from  orbit  to  orbit.  Since  in  these  formulas  the  perturbing  influence  of 
the  oblateness  of  the  planet  is  not  taken  into  account,  the  computed  results  will 
be  close  to  the  real  picture  of  the  evolution  of  very  high  orbits  with  uniform 
physical  characteristics  of  the  satellite  and  for  lighter  satellites  during  single 
orbits. 

In  this  paper  we  investigate  the  Long-period  perturbations  of  the  satellite 

orbital  elements.  We  shall  not  consider  short-period  perturbations  caused  by 

shadowing  effects.  The  sequence  for  comparing  the  evolution  of  the  orbits  is  as 

2 

follows:  we  first  set  up  the  model  of  a  light  spherical  satellite  (A/m  *  200  cm  /g) 
and  consider  the  perturbations  for  orbits  of  different  dimensions  and  orientations 
around  the  earth;  second,  we  consider  different  satellites  to  study  the  influence 
of  the  parameter  A/m;  third,  we  examine  two  identical  satellite  orbits  for  Earth 
and  Mara  (with  like  parameters  A/m)  to  compare  the  satellite  perturbations  for  differ¬ 
ent  planets. 

-4  2 

We  now  turn  to  perturbations  of  earth  satellites  (P  **  0.46  x  10  g/cm*sec  ) 
and  examine  aeparctely  the  perturbations  of  various  periods.  It  follows  from  the 
form  of  formula  (4)  that  the  perturbation  periods  are  determined  by  the  rates  of 
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change  of  the  quantities  K^,  jl:,  i,  0.  Obviously,  «e  must  first  of  all  consider 
the  yearly  periodicity  related  to  solar  motion;  this  is  conveniently  done  by  vary¬ 
ing  the  elements  e  and  u),  which  are  subjected  to  the  largest  relative  changes  due 
to  radiation  pressure  of  the  elements.  The  elements  i  and  Q  are  changed  rather 
weakly  over  the  course  of  a  year  (especially  with  small  values  of  e) ,  since  e 
appears  linearly  in  the  coefficient  in  the  year  term  in  Q^.  The  semi-major  axis  a, 
in  general,  does  not  change  with  full  illumination.  The  change  in  is  determined 
basically  by  the  change  in  Ji  and  will  not  be  considered  individually. 

Dependence  on  Semi-major  Axes.  Let  us  consider  an  orbit  vith  the  elements 
e0  3  0.1,  L,  *  5° ,  0-  *  JCj  =  3  0°.  Figs.  1  and  2  show  graphically  the  changes  in 

e  and  &  over  a  year  for  the  following  values  of  a:  24,000,  21,400,  13,500,  7.400  km, 
as  computed  using  formula  (5).  The  curves  of  e(a)  shew  the  annual  oscillation,  with 
the  curve  changing  from  convex-upward  to  convex-downward  (the  results  of  the  calcu¬ 
lations  for  a  3  42,000  km  are  in  agreement  with  the  data  of  references  [4j  and  L 5 j ) . 
As  the  computations  for  other  initial  data  showed,  as  eg  increases  the  entire  set  of 
curves  will  have  a  tendency  to  fall  lower.  For  each  value  of  ec  it  is  possible  to 
select  a  value  of  a  fox  which  a  minimum  change  in  e  is  observed  (for  example,  0.1  * 
e  £  0.2,  13,000  km  s  a  £  42,000  kra) .  The  shape  of  the  eccentricity  curves  repeats 
from  year  to  year.  The  perigee  (Fig. 2)  moves  in  an  rnnual  circle  with  variable 
velocity,  making  a  full  revolution  per  year  (the  straight  line  denotes  the  change 
in  longitude  of  the  sun).  With  large  values  of  eg  or  small  values  of  A/m  the  per¬ 
turbations  in  us,  as  follows  from  (5),  will  be  small  and  will  show  up  as  oscillatory 
(loop-shaped)  motion  around  the  initial  position. 

Dependence  on  Eccentricity.  Let  us  select  an  orbit  close  to  synchronous: 
a  *  42,000  km,  i,*  5°,  ui,  *  O,  =  \%  *  o'5.  Figs.  3  and  4  show  graphs  of  e  and  u>  for 
a  series  of  e„ :  0.01,  0.1,  0.2,  0.3,  0.4  and  0.5.  Again  we  see  three  types  of  curve 
for  e(e0):  for  e0  s  0.32  the  perigee  moves  in  a  forward  motion;  for  e  >  0*32  the 
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notion  reverses.  For  «„  «  0.5  the  perigee  moves  with  a  counter  loop-shaped  motion, 
making  k  complete  revolution  in  22  years.  Fig.  5  illustrates  such  s  perigee  path 
over  a  period  of  4  years:  points  A^  and  are  separated  by  year  intervals,  A^  is 
the  sii?'  1  position  at  these  times.  As  will  be  shewn  below,  the  perigee  motion  for 
the  next  22  years  will  be  exactly  the  opposite  since  the  total  period  of  change  is 
44  years  (*or  a  synchronous  satellite). 

Influence  Inclination.  For  orbits  with  a  *  42,000  km.  a-,  *  0.01, 
fij  =  v'r  *  0°  we  examine  a  series  of  inclinations  i0  *  1° ,  23c ,  46°,  70^,  90° , 

11  ,  ’SG°,  160°  and  17$°.  AH  curves  of  e(iu)  (Fig. 6)  are  distributed  In  a  region 

..  *  *i<  by  the  two  curves  £<1.,  a,  c  «  23°)  and  e(i2  *=  90°  +  t  «  113°),  where  the 

for  values  of  i0  symmetrical  with  respect  to  i^  or  must  coincide.  This 
regularity  is  observed  for  any  choice  of  initial  elements,  although  the  shape  of 
the  “(i^)  curves  will  in  this  case  be  altered. 

Effect  of  Nodal  Longitude  a-d  Argument  of  Perigee.  Let  us  select  two  iden¬ 
tical  orbits:  a  *  42,000  km,  *»  O.i,  ^  *  23.5°,  >£.  -  0°.  The  effects  of  the  annual 
variation  in  i  and  Q  is  completely  excluded  from  the  perturbations  of  e  and  x  for 
such  an  orbit,  since  for  %  *  0? ,  -51  «  $0  *  OP  (ecliptic  orbit),  and  these  per¬ 
turbations  are  small  for  £$,  f  & .  We  shall  choose  the  values  =  0s ,  (1,  »  & ,  90°, 

1SCP,  270°  (Fig.  7)  and  0,  -  OP,  0°,  9(P,  18<f,  27 OP  (Fig. 8).  It  is  apparent 

that  the  evolution  of  the  elements  is  not  affectea  by  the  actual  values  of  Qj  or 
JX,,  but  by  the  mutual  initial  disposition  of  node  and  perigee.  The  inequality  of 
amplitudes  in  Figs.  ?  «nd  8  indicates  that  Q  and  u»  enter  in  quite  different  ways 
into  the  equation. 

Influence  of  Solar  Longitude.  Let  us  consider  a  synchronous  orbit:  a  * 

42,000  km,  e0  ■*  0.1,  tj  ~  23.  S*5  (eclipcie),  ■  %  "0s.  The  curves  of  e  and  w 

(\®  *  G,  18CP)  coincide,  and  the  curves  of  e  and  (\%  *  90°,  270^ }  are  symmetri- 
caily  disposed  relative  to  the  straight  line  t  *  P/2  (Figa.  9  and  10).  In  the 
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general  case  of  arbitrary  eccentricity  the  curves  for  \°  =  3°  and  180°  are  symmetric 
cally  disposed  relative  to  the  straight  line  e0  =  const.  On  the  whole  the  asymmetry 
is  determined  by  the  asymmetry  in  the  positions  of  perigee  and  the  sun. 

Let  us  now  examine  the  longer-period  perturbations.  Here  we  shall  consider 
the  behavior  of  the  elements  i  and  ft  when  the  amplitudes  of  the  corresponding  varia¬ 
tions  of  e  and  a s  are  small  compared  to  their  annual  amplitude  (the  ordinates  of 
the  maxima  of  e  exhibit  only  a  weak  long-period  trend).  The  annual  changes  in  e, 
jj,  and  \  exert  an  influence  on  i  and  Q  in  the  form  of  weak  annual  variations 

v> 

superposed  on  the  main  trend.  These  variations  will  not  be  shown  on  the  following 
graphs . 


The  variation  of  i  and  Q  without  taking  shadowing  info  account  is  defined 
by  the  formulas 


<U_ 

At 


3  ncrci 

Y/7^7i 


Wcos«, 


AQ  3  na~n 

M  2  sin 


W  sin  w, 


(7) 


where 

Wz=  sin  /  sin*  sin  (2  -f-  X0)  -j-  sin  i  cos*  4  sin  (2  —  X3)  -j- 
+  cos  i  sin  s  sin  X©. 

It  is  clear  from  (7)  that  we  have,  as  a  function  of  the  arguments  of  the  trigono¬ 
metric  terms,  either  the  annual  variations  due  to  changes  in  or  the  variations 

tii.  clQ 

arising  from  the  superposition  of  variations  with  periods  2n:  —  and  2*t:  ~  .  Since 
these  derivatives  are  small  the  oscillation  periods  can  be  quite  long. 
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10  20  30  hG  50  60  70 
Fig. 11 


10  20  30  40  50  CO  70 
Fig. 12 


Influence  of  Seml^axis.  We  choose  an  elliptical  orbit,  near  equatorial 
(e0=  0.1,  1q  “  5°  ,  Jtf0  ■  (\  =  =  OP)  and  vary  the  semiaxis:  a^  =  42,000  km,  * 

23,000  km,  a^  *  13,500  km.  From  the  formulas  it  follows  that  the  amplitudes  should 
be  related  as  /«"  and  the  periods  as  1/  /a. 

Graphs  of  the  variations  in  i  and  C.  are  shown  in  Figs.  11  and  12.  The 
amplitudes  are  related  as  2.5:2: 1.5,  the  periods  as  4:5. 2:6. 7.  The  minimum  period 
represents  a  stationary  orbit  --  44  years.  The  variations  in  the  other  elements  will 
be  considered  just  for  this  period. 

The  Influence  of  Elements  having  an  Annual  Period  of  Variation.  Among 
these  quantities  are  e,  x  and  X^.  By  virtue  of  the  fact  that  has  no  annual  per* 
iodicity  and  it  is  negligibly  small  for  e  and  i>,  these  elements  exert  no  appreciable 
effect,  on  the  course  of  i  and  w.  Calculations  for  the  stable  orbit  yield  almost 
identical  curves  and  therefore  we  shall  not  reproduce  them  here.  This  is  observed 
even  for  large  values  of  e,  when  the  annual  movement  in  ju  is  appreciable. 


»  20  30  W  50 

t,  years 


Fig. 13 


10  20  33  « 

t,  years 

Fig. 14 
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Influence  of  Inclination.  The  elements  i  and  0  react  strongly  on  the 
variation  in  L,  and  their  plots  differ  for  different  values  of  ^ .  Let  us  again 
consider  the  orbit;  a  *  42,000  km,  eQ  ■  0.01,  •*  “  A^  ■  0° .  All  equatorial 

and  near-  equatorial  orbits  show  an  increase  in  i  proportional  to  1/sin  i  (Fig.  13) 
until,  finally,  the  curve  of  i  (i^  as  e)  reduces  to  a  straight  line  (1J  *  0, 

&i  *  50  -  0).  For  all  ic  >  e,  i  begins  to  decrease,  with  the  maximum  amplitude  of 

the  decrease  corresponding  to  as  9 (P  +  e  «  i!3° ,  since  W(1q  as  90°  +  e)  * 

cos  £  sin  A  (cos  e  +  sin  e)  has  a  maximum.  Then  the  amplitudes  decrease  and  the 
curves  i(ig  *  1?°°)  and  1(1^  •»  1°)  are  symmetrical  around  the  line  *  90°.  The 

curves  of  0(io)  (Fig.  14)  show  a  long-period  counter  motion  of  the  node  for  >  e 

and  a  loop-shaped  motion  for  <  e .  The  node  of  the  elliptical  orbit  is  stationary. 
The  distribution  of  curves  by  amplitude  is  preserved  for  any  value  of  but  the 
symmetry  is  lost. 

Influence  of  Node  Longitude.  Let  us  select  a  stationary,  ecliptic  orbit. 
In  this  case  we  observe  complete  symmetry  of  the  curves  i  (Q,  *  0°,  180°)  for  both 
half-periods,  where  5i(0  »  0)  ■  0  since  W  ■  0.  The  curves  i(Qj  *  90°,  270°)  pass 
through  each  other  on  rotation  around  the  axis  t  *  P/2  (Fig.  15).  In  the  case 
«*  0°  we  again  have  50  *  0  (Fig.  16),  and  for  other  values  of  %  the  node  motion 
is  in  the  opposite  direction  with  variable  velocities.  The  entire  set  of  curves 
i(C^)  for  equatorial  orbita  should  be  stretched  out  upward,  and  the  rate  of  change 
of  0(03 )  should  vary  nsore  strongly.  In  the  case  of  a  polar  orbit  the  curves 
i(Qj  ■  0°,  180°)  should  be  rigorously  symmetrical  with  respect  to  Lj  *>  const,  end 
the  rate  of  movement  of  the  node  should  be  almost  constant. 
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Let  us  now  examine  the  influence  of  the  parameter  A/m  on  the  development 

of  the  orbital  elements.  Until  r.ow  we  have  varied  the  orbital  elements  of  earth 

2 

satellites  vith  the  same  value  of  A/m  =  200  cro'/g,  a»id  have  considered  this  to  be 
a  fixed  quantity.  Let  us  now  set  up  a  single  orbit  (a  =  42,000  km,  e5  =  0.01, 

L,**  5° ,  %  •  Qj  ■  0°)  and  select  certain  values  of  A/m  (A/ra  =  200:0.1  A/m, 

5  A/ra,  7  A/ra,  10  A/ra).  Fig.  17  shows  curves  of  the  change  in  eccentricity  over  a 


Fig. 18 

2 -month  interval.  From  these  curves  we  can  form  a  conclusion  regarding  the  magni¬ 
tude  of  the  changes  in  the  remaining  orbital  elements.  In  the  calculations  we  used, 
as  before,  the  spherical  satellite  model  with  constant  A/m.  For  this  satellite 
the  question  of  orientation  relative  to  the  solar  radiation  during  the  course  of 
the  year  does  not  arise.  In  addition,  the  acceleration  of  a  spherical  satellite 
does  not  depend  on  its  surface  reflectivity  (with  geometrical  reflection  properties). 
In  the  general  case  the  object  with  constant  A/m  may  be  a  satellite  of  arbitrary 
shape  that  preserves  its  orientation  relative  to  the  sun.  If  the  orientation  of 
the  satellite  is  ".rbitrary,  the  law  of  the  change  in  illuminated  surface  should  be 
given  functionally,  33  was  done,  for  example,  in  reference  [5j  for  synchronous  sate¬ 
llites.  Here  we  must  keep  in  mind  that  the  acceleration  F  of  any  satellite  wj.ll  lie 
in  the  range 

Fg  <F<Fp 


•here  V  la  the  acceleration  of  a  sphere  and  F^  is  the  acceleration  of  a  flat  plate 
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Hi 


whOi'e  illuminated  surface  is  equal  to  the  central  cross  section  of  the  sphere  and 
of  the.  arbitrary  satellite.  In  this  case  we  compute  the  acceleration  from  the 


forma'.  »s 

=p,4.  =vn+R. 

Z7.  — (1  -f*c)P.cos,d“1  Fy  =(1  —  s)  P,  s!n  t)  C05  0— , 

where  £  is  the  coefficient  of  reflection,  is  the  angle  of  incidence,  and  FT 

are  the  normal  and  tangential  components  of  the  pressure. 

In  order  to  compare  the  perturbations  of  individual  satellites  for  two 

planets,  let  us  consider  satellites  orbiting  Earth  and  Mars.  For  comparison  pur- 

2 

poses  we  select  a  heavy  satellite  (A/a  *  0.1  cm  /g)  in  an  orbit  with  the  elements 
a  =  42,000  km,  e,  «  0.01,  i0  =  5°  ,  «%=.%  *  X°  ~  0°.  For  the  earth  this  is  close 
to  an  actual  synchronous  satellite.  Fig.  18  shows  curves  of  the  change  in  the 
perigee  radius-vector  in  the  course  of  a  year.  The  maximum  amplitude  for  the  earth 
amounts  to  about  7  km.  The  corresponding  amplitudes  in  uj  are  ml0  per  year;  the 
perigee  executes  small  oscillations  around  the  stable  position,  so  that  a  longitude 
shift'  of  the  satellite  is  not  expected.  The  middle  curve  of  Fig.  18  shows  for 
a  satellite  of  Mars;  the  annual  amplitude  does  not  exceed  0.5  km.  This  difference 
is  explained  by  the  fact  that  the  value  of  Pr  in  the  Mars  orbit  is  about  half  that 
in  the  Earth  orbit,  which  is  close  to  the  synchronous  orbit  (a  s.-  20,000  km).  We 
present  a  table  of  some  values  taken  from  the  calculations  made  for  the  curves  of 
Fig.  18. 


Earth 

Mars 

'•  cra.sec£ 

i 

0.45-  JO-4 

0J20-10— 4 

F.cm/sec^ 

L  * 

0.4S-!0~ * 

0.20- i«r5 

R,  KM 

6371 

3407 

fM,  Ki&i  sec' 

358600 

42900 

«• 

23c5 

25'0 

iQ,  deg/d*y 

0.S8565 

0.52500 
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